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Several teams fiave reported peculiar frequency spectra for flows in a spherical shell. To 
address their origin, we perform numerical simulations of the spherical Couette flow in 
a dipolar magnetic field, in the configuration of the DTS experiment. The frequency 
spectra computed from time-series of the induced magnetic field display similar bumpy 
spectra, where each bump corresponds to a given azimuthal mode number m. The bumps 
show up at moderate Reynolds number {c^ 2 600) if the time-series are long enough 
(> 300 rotations of the inner sphere). We present a new method that permits to retrieve 
the dominant frequencies for individual mode numbers m, and to extract the modal 
structure of the full non-linear flow. The maps of the energy of the fluctuations and 
the spatio-temporal evolution of the velocity field suggest that fluctuations originate in 
the outer boundary layer. The threshold of instability if found at Rcc = 1 860. The 
fluctuations result from two coupled instabilities: high latitude Bodewadt-type boundary 
layer instability, and secondary non-axisymmetric instability of a centripetal jet forming 
at the equator of the outer sphere. We explore the variation of the magnetic and kinetic 
energies with the input parameters, and show that a modified Elsasser number controls 
their evolution. We can thus compare with experimental determinations of these energies 
and find a good agreement. Because of the dipolar nature of the imposed magnetic field, 
the energy of magnetic fluctuations is much larger near the inner sphere, but their origin 
lies in velocity fluctuations that initiate in the outer boundary layer. 
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1. Introduction 

It is now well established that the magnetic field of mo st planets and stars is generated 
by the dynamo mechanism ( Larmorlll919l : [Elsasser 1946[ ). Motions within an electrically 
conducting medium can amplify infinitesimally small magnetic field fluctuations up to 
a level where the Lorentz force that results is large enough to stop their amplification. 
This is possible for large enough values of the magnetic Reynolds number Rm = UL/r] 
(where is a typical flow velocity, L a typical length, and r] is the magnetic diffusivity 
of the mediun i) . 

Analytical ( BusselfToyoh and numerical ( Glatzmaier fc Robertsl[l995 ) convective dy- 
namo models, in which the flow is driven by the buoyancy force of thermal or composi- 
tional origin, have demonstrated the relevance of the dynamo mechanism for generating 
the Earth's magnetic field. Other forcings, due to prece ssion, tides or impacts are also 
invoked to explain the fields of some other planets (Le Bars et al\\2011\ ). 
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In year 2000, two experime nts demonstrated dynamo action in the lab (jGailitis et 



2OOII : IStieglitz fc Mullerll200ll ). In both cases, the forcing was mechanical, with a domi- 



nant large-scale fl ow. Efforts to produce dynamo action with a highly turbulent fl ow are 
stih on the way (|LathroD fc ForestI I2OI1I : iKaplan et alll201ll : iFrick et a/Jl2010[ ). while 
a rich variety of dynamo behaviours have been discovered in the VKS experiment 



(jBerhanu et a/.l 120071 : iMonchaux et a/.l 120071 ) when ferromagnetic disks stir the fluid. 

All these experiments use liquid sodium as a working fluid. The magnetic Prandtl 
number Pm = u/r] of liquid sodium is less than 10~^ {u is the kinematic viscosity), 
so that experiments that achieve Rm of order 50 (as required for dynamo action) have 
kinetic Reynolds number Re = UL/v in excess of 10^. This contrasts with numerical 
simulations, which require heavy computations with 1024^ grid points to reach Re = 
10^. Since Reynolds numbers of flows in planetary cores and stars are much larger, we 
have to rely on theory to bridge the gap. Dynamo turbulence is a crucial issue because 
dissipation is very much dependent upon the scale and strength of turbulent fluctuations. 
The question of instabili ties and turbulence is also central in the study of accretion disks 
( Balbus Hawlevlll99ll) . Laboratory experiments can bring some constraints since they 
exhibit intermediate Reynolds numbers. 

In that respect, the observation in several experiments of very peculiar frequency spec- 
tra, characterized by a succession of peaks or bumps deserves some attention. Such bumpy 
spectra have been obtained in both spherical and cylindrical geometries, when rotation 
or/and magnetic flelds are present, two ingredients that also play a major role in natural 
s ystems. 

iKellev et all (|2007[ ) were the flrst to observe a bumpy spectrum in a rotating spherical 
Couette experiment. A small axial magnetic fleld was applied and the induced fleld 
was used as a marker of the flow. The authors showed that the frequency and pattern 
of the modes correspond to a set of inertial modes. Inertial modes are the oscillatory 
linear response of a fluid to a time-dependent perturbation where the Coriolis force is 
the restoring force. Several hypotheses have been put forward to explain the excitation 
of ine rtial modes in the se experiments: over critical reflection off the inner Stewart son 
l aver (|Kellev et 



2010[ ). and turbulence from the tangent cylinder on the inner sphere 



Matsui et a/.ll201l[ ). 

Bumpy frequency spectra were also report ed by Schmitt a/1 (l2008h in the j^T^S mag- 
netized spherical Couette flow experiment ( Cardin all 200a iNataf et a/.l I2OO6U2OO8I : 
Brito et a/1l201l[ ). An example is shown in flgure[lK. Schmitt et all ( 2008[ ) could show, by 



correlating signals measured at several lon gitudes, that each bu mp is characterized by a 
given azimuthal wavenumber m (flgure[Tl3). Schmitt et all ( 20121 ) further investigated the 
properties of the bumps and showed a good correspondence with linear magneto-inertial 
modes, in which both the Coriolis and the Lorentz forces play a leading role. 

Finally, modes of azimuthal wavenumber m = 1 were observed in two magnetized 
Couette flow experiments aimed at detecting the magneto-rotational instability (MRI): 
in spherica l geometry in Maryla nd (Sisan et al. 2004) and in cylindrical geometry in 
Princeton 



isan et all (|2004l ) discovered magnetic modes that 



v laryla nd ( 
I l2010[ ). [S^ 



appeared only when the imposed magnetic fleld was strong enough and interpreted their 
observations as evidence for the MRI, even though the most unstable mode is expected 
to be axisymmetric (m = 0) in their geometry. Rotating spheric al Couette flow in an 
axial magnetic fleld was studied numerically by iHollerbcichl (|2009[ ) , who suggested that 
instabilities of the meridi onal circulat i on in t he equatorial region co uld account for some 



of the modes observed by Sisan et~all (|2004{ ). ICissinger et al. ( 2011 ) further investigated 



this situation and showed that the instabilities that affect the Stewartson layer around 
the inner sphere, modifled by the imposed magnetic fleld, have properties similar to the 
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Figure 1. Spectral bumps in the DTS magnetized spherical Couette flow experiment, a) Fre- 
quency spectra of the three components of the surface magnetic field (5r, be, b^ in this order 
going up) measured at two different latitudes (5° and 45° as indicated). The time window used 
lasts 4 000 turns. Frequencies (x-axis) are normalized by the rotation frequency of the inner 
sphere (/ = —13.3 Hz). The spectral energy density (^-axis) is normalized as in section [4Jl to 
facilitate the comparison with the corresponding simulation results. The vertical scale applies 
to the lowest spectrum and successive spectra are shifted by one decade for clarity (two decades 
between different latitudes). Note the succession of bumps that dominate the spectra, b) co- 
variance between two br time-series recorded at points 128° apart in longitude (same latitude 
= —35°), in a frequency- azimuthal mode number (m) plot. This plot shows that each spectral 
bumps in a) corresponds to a well-defined integer m, which increases with frequency (step- wise 
succession of positive covariance values starting at m = 0). 



MR L In contrast to the st andard MRI, the instabilities evidenced bv lHollerbachl (|2Q09[ ) 
and iGissinger et all ( 2Qlll ) are induct ionless. 



In the cylindrical Taylor- Couette geometry with a strong imposed axial field, (iNornberg et al 
2QlQh observed m = 1 rotating modes. They claimed that these modes could be iden- 



tified to the fast and slow magneto- Cor iolis waves expected to develop when both the 
Coriolis and Lorentz forces have a comparable strength. Considering the fast magnetic 
diffusion in their experiment ( Lundquist number of about 2), this interpretation was 



rather questionable, and indeed [Roach et 'al. ( 2012 ) have recently reinterpreted these ob- 



se rvations in term s of in stabilities of an internal shear layer, in the spirit of the findings 
of IGissinger et al\ (|201l[ ). 

Clearly, magnetized Couette flows display a rich palette of modes and instabilities, 
and it is important to identify the proper mechanisms in order to extrapolate to natural 
systems. Hollerbach has investigated the instabilities of magnetized spherical Couette flow 
in a series of numerical simulations ( Hollerbach fc Skinner! 12001; Hollerbach et al. ^ l2007[ 
Hollerbachl[2009l ). However, bumpy spectra as observed by ^chmitt et al. (2008) were 



never mentioned. In this article, we perform numerical simulations in the geometry of 
the DTS expe riment, and focus on the origin of these bumpy spectra. The observations of 
Schmitt et al\ (j2008) are illustrated by figure [H but the reader should refer to their article 



for a more detailed presentation. More specifically, we wish to answer the following major 
questions: how and where are the various modes excited? Are these spectra observed 
because of the large value of the Reynolds number? 
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Figure 2. Sketch of the DTS experiment as modeled in this article. The inner sphere of radius 
ri rotates around the vertical axis at angular velocity Q. It consists in a copper shell enclosing 
a permanent magnet, which produces the imposed dipolar magnetic field Bd- Liquid sodium of 
electric conductivity a fills the gap between the inner sphere and the stainless steel outer shell 
of inner radius To- The spherical coordinate system we use is drawn. 

We present the numerical model and the mean flow in section [2j We perform spec- 
tral analyses in section [3l and investigate fluctuations and instabilities in section HI A 
discussion concludes the article. 

2. Numerical model and mean flow 

The DTS experiment that we wish to model is a spherical Couette flow experiment 
with an imposed dipolar magnetic field. Liquid sodium is used as a working fluid. It is 
contained between an inner sphere and a concentric outer shell, from radius r = to 
r = To {vi = 74 mm, Tq = 210 mm). The inner sphere consists of a 15 mm-thick copper 
shell, which encloses a permanent magnet that produces the imposed magnetic field, 
whose intensity reaches 175 mT at the equator of the inner sphere. The stainless steel 
outer shell is 5 mm thick. The inner sphere can rotate around the vertical axis (which is 
the axis of the dipole) at rotation rates / = 27rl] up to 30 Hz. Although the outer shell 
can also rotate independently around the vertical axis in DTS, we only consider here the 
case when the outer sphere is at rest. 

All these elements are taken into account in the numerical model, which is sketched 
in figure [21 In particular, we reproduce the ratio in electric conductivity of the three 
materials (copper, sodium, stainless steel). In the experiment, the inner sphere is held by 
25 mm-diameter stainless steel shafts, which are not included in the numerical model. 

2.1. Equations 

We solve the Navier- Stokes and magnetic induction equations that govern the evolution 
of the velocity and magnetic fields of an incompressible fluid in a spherical shell: 




-Vp + z/V^u H (V X B) X B, 

P Pop 



(2.1) 



— - = Vx{uxB)-Vx {r]{r)V x B) , 



(2.2) 
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V.n = 0, V'B = 0, (2.3) 

where u and p stand for the velocity and pressure fields respectively. Time is denoted 
by while p and u are the density and kinematic viscosity of the fluid. The magnetic 
diffusivity r]{r) is given by r]{r) = {fiQa{r))~^ where a{r) is the electric conductivity of 
the medium (fluid or solid shells) and /io the magnetic permeability of vacuum. In the 
fluid, the conductivity cr(r) = cr^a is constant. The last term of the flrst equation is the 
Lorentz force. B is the magnetic fleld. It contains the imposed dipolar magnetic fleld Bd 
given by: 

3 

Bd = Bo^ (2cosi9e^ + sin6>e0) , (2.4) 

where is the colatitude, Cr and bq are the unitary vectors in the radial and orthoradial 
directions. Bq is the intensity of the fleld at the equator on the outer surface of the fluid 
(r = To). 

2.2. Boundary conditions 
We use no-slip boundary conditions for the velocity fleld on the inner and outer surfaces: 

u = QrsmOe(p for r ^ r^, u = for r ^ ro- (2-5) 

We model the copper shell that holds the magnet in DTS as a conductive shell with 
electric conductivity gcu = ^'"^cFNa- The outer stainless steel shell is modelled as a shell 
of conductivity gss = cFNa/^- These values reproduce the experimental conductivity 
contrasts. The conductivity jumps are implemented by taking a continuous radial con- 
ductivity proflle with sharp localized variations at both interfaces (3 to 5 densifled grid 
points). The internal magnet and the medium beyond the outer stainless steel shell are 
modelled as electric insulators. The magnetic fleld thus matches potential flelds at the 
inner and outer surfaces. 

2.3. Numerical scheme 

Our three-dimensional spherical code (XSHELLS) uses second-order flnite differences in 
radius and pseudo-spectral spherical harmonic expansion, f or which it relie s on the very 
efficient spherical harmonic transform of the SHTns library (ISchaefferll2QT2l ). It performs 



the time-stepping of the momentum equation in the ffuid spherical shell, and the time- 
stepping of the induction equation both in the conducting walls and in the ffuid. It uses a 
semi- implicit Crank-Nicholson scheme for the diffusive terms, while the non-linear terms 
are handled by an Adams-Bashforth scheme (second order in time). The simulations that 
we present typically have 600 radial grid points (with a signiffcant concentration near 
the interfaces) while the spherical harmonic expansion is carried up to degree 120 and 
order 40. 

2.4. Dimensionless parameters 

We deflne in table [1] the dimensionless numbers that govern the solutions in our problem. 
We pick the outer radius as a length scale, and the intensity of the magnetic fleld 
at the equator of the outer surface, as a magnetic fleld scale. Note that, due to the dipolar 
nature of the imposed magnetic fleld, its intensity is 23 times larger at the equator of 
the inner sphere. The angular velocity of the inner sphere yields the inverse of the time 
scale. We choose U = l^r^, the tangential velocity at the equator of the inner sphere, as 
typical velocity. 
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symbol expression simulations simulations DTS Earth's 

fixed Ha fixed Re / = [3 — 30] Hz core 

Pm ly/r] 10'^ 10"^ 7 x 10"^ 10"^ 

Re nnro/iy [2 611 - 10 100] 2 611 [0.45 - 4.5] x 10^ 10^ 

A aBl/pVt [3.4 - 0.9] x 10"^ [0.21 - 13.1] x 10"^ [3 - 0.3] x 10"^ 10 

Ha roBo/y/pfioi^r] 16 [4-31] 200 10^ 

Table 1. Typical values of the dimensionless numbers in the numerical simulations, in the DTS 
experiment (computed for / = Q/27r = 3 — 30 Hz) and in the Earth's core. The dimensionless 
numbers of our reference simulation are in bold. 



The solutions are governed by three independent dimensionless numbers but several 
combinations are possible and we try to pick the most relevant ones. The magnetic 
Prandtl number Pm compares the diffusion of momentum to that of the magnetic field. 
It is small in both the simulations and the experiment. The Reynolds number Re is 
of course essential, as it determines the level of fluctuations. It is not feasible to run 
numerical simulations with Reynolds numbers as large a s in the DTS experiment. 

However, one of the main findings of lBrito et all (|201ll ) is that, because of the imposed 



dipolar magnetic field, the time-averaged flow is mainly governed by the balance between 
the Lorentz and the Coriolis forces, where the latter is due to the global rotation of 
the fluid, which is very efficiently entrained by the inner sphere, ev en when the outer 



sphere is at rest. That balance is measured by the Elsasser number A. lBrito et al\ (|201lf ) 
showed that one can recover the proper balance at achievable values of Re by reducing 
the influence of the magnetic field, keeping the effective Elsasser number A as in the 
experiment. 



Neverthe less, Cardin et all ( 20021 ) introduced another number A (named the Lehnert 



number bv ljaultl (|2008l )). which provides a better measure of this balance for fast time- 
dependent phenomena. The Lehnert number A compares the periods of Alfven waves to 
that of inert ial waves. It is given by: 

^=W^- (2.6) 

In the Earth's core, this number is of order 10~^ and inertial waves dominate. They 
force the flow to be quasi-geostrophic on short time-scales (|jaultl[2QQ8l ). However, mag- 
netic diffusion severely limits the propagation of Alfven waves in the DTS experiment. 
This is measured by the Lundquist number, which is the ratio of the magnetic diffusion 
time to the typical transit time of an Alfven wave across the sphere, here given by: 

Lu = — ; , (2.7) 

which is taken as Lu = 0.5 in the numerical simulations, in agreement with the exper- 
imental value. 



We therefore follow the same strategy as iBrito et al\ (|201ll ) , and try to keep the El- 



sasser number of the numerical simulations similar to its experimental value. Our refer- 
ence case thus has: Pm = 10~^, Re = 2 611 and A = 3.4 x 10~^. The Hartmann number 
is Ha = 16, quite smaller that its experimental counterpart {Ha = 200). It follows that 
A = 6.8 X 10~^ and Rm = RePm = 2.6 for the reference case. Most results shown in this 
article relate to our reference case, but we also present some results computed for other 
Reynolds and Hartmann numbers, as indicated in tabled) 



Modes and instabilities in magnetized spherical Couette flow 7 




Figure 3. Time- averaged meridional slice of the velocity field, a) Angular velocity isovalues. 
Note the zone of super-rotation near the inner sphere, b) Meridional streamlines. The maximum 
meridional velocity is 0.4. Pm = 10"^, Re = 2 611, A = 3.4 x 10"^ 



2.5. Mean flow 

The time-averaged properties of the m agnetized spherical Couette flow have been inves- 
tigated in detail bv lBrito et al. (|2Qll[ ). and we simply recall here a few key observations. 
We plot in figure [3] the time- averaged velocity field in a meridional plane, for our reference 
simulation (Pm = 10~^, Re = 2 611 and A = 3.4 x 10~^). Two distinct regions show 
up in the map of mean angular velocity (figure [3^): an outer geostrophic region, where 
the angular velocity only v aries with the cylindrical radius s; an inner region that obeys 
Ferraro law (|Ferrardll937[ ) around the equator: the angular velocity is constant along 
field lines of the imposed dipolar magnetic field. Note the presence of a thin boundary 
layer at the outer surface. The poloidal streamlines (figure [Sb) display a circulation from 
the equator towards the poles beneath the outer surface, where the polewards velocity 
reaches 0.4 Qri. 

To check our numerical set-up, we compare the time- averaged velocity field of our 
simulation with that obtained by iBrito et all ( 201ll ) using an independent axisymmet- 
ric equatorially-symmetric code. The parameters and boundary conditions are identical, 
except that th e magnetic bounda ry condition at r = is treated in the thin-shell ap- 
proximation in lBrito et al\ ( 201ll ). 

Figure |4] c ompares the radial profile of the angular velocity in the equatorial plane 
computed by iBrito et o/l ( 201l[ ) to our axisymmetric solution, averaged over 50 rota- 
tion times, and to our 3D spherical solution, averaged over 100 rotation times. The 
two axisymmetric solutions agree very well, while the 3D solution exhibits a slightly 
lower angular velocity near the outer surface. Note that the angular velocity of the fluid 
reaches values as high as 20% larger than that of t he inn er sphere. This phenomenon of 
super-rotation wa s first predicted bv lDormv et al\ ( 1998 ) in the same geometry (also see 
Starchenkol ( 19971 )). but in their linear study, the zone of super- rotation was encl o sed in 
the magnetic field line touching the equator of the outer sphere. Hollerbach et all (|2007l ) 
showed that non-linear terms shift the zone of super-rotation from the outer sphere to 
close to the inner sphere, as observed here. The excess of 20% is in good agreement with 
the super-rotation measured in the DTS experiment for / = 3 Hz {Re = 4.5 x 10^) 



(jBrito et a/Jl201lf ). 
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Figure 4. Radial profiles of time- averaged angular velocity in the equatorial plane for different 
simulations with the same parameters (Pm = 10"^ , Re 2611 and A = 3.4 x 10"^). Dotted line: 
axisymmetric equatorially-symmetric solution of iBrito et aL\ (|2011l l: solid line: axisymmetric 
solution computed with our XSHELLS code; dashed line: 3D solution from XSHELLS. The 
curves at the bottom give the unsigned difference between the 3D solution and the axisymmetric 
one, and between the two axisymmetric solutions, scaled up by a factor 5. 
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Figure 5. Spectral bumps in our reference numerical simulation {Pm 
and A = 3.4 x 10~^). a) Frequency spectra of the three components of the magnetic field (br, 
he, in this order going up) recorded at two different latitudes (0° and 41° as indicated). 
Same adimensionalization and plotting conventions as in figure [TJi. The time window used lasts 
540 turns. Note the spectral bumps and compare with figure [TJi. b) m versus frequency plot 
for the same run, obtained from the covariance of two time-series recorded at points 72° 
apart in longitude (same latitude = 41°). The red patches for negative m indicate that the 
successive spectral bumps have a well-defined mode number m, whose absolute value increases 
with frequency (compare with figure [TJ) where the relevant m are positive because the inner 
sphere spins in the negative direction). 



3. Spectra and modes 

3.1. Frequency spectra 

In order to compare the numerical results to the experimental measurements of the 
fluctuations, we perform a simulation over a long time- window (600 rotation periods) and 
record the magnetic field induced at the surface at selected latitudes. We then compute 
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^shifted 



Figure 6. Partial energy frequency spectra Smi"^) of the equat or ially- symmetric magnetic field 
at r < 0.55. Frequencies are normalized by /, the rotation frequency of the inner sphere (vertical 
dashed line), a) Raw spectra for azimuthal mode numbers m = to 4 (shifted vertically for 
clarity). The vertical solid line indicates the frequency ui* at which we retrieve the m — 2 
mode structure in figure [51 b) Same as a) except that frequencies are shifted according to: 
^shifted = taj — mf fluid (see text). 



the power spectra of these records as a function of frequency. Typical spectra are shown 
in figure [5^. A sequence of bumps is clearly visible for both the radial and the azimuthal 
components of the magnetic field. The spectra do not display power-law behaviour. 

We note that long time series (longer than 300 rotation periods) are needed for the 
spectral bumps to show up clearly. The bumps are not as pronounced as in figur^Ik, but 
we note that 4 000 turns were used for these experimental spectra. It could also be that 
the bumps are enhanced at higher Reynolds number. 



3.2. Azimuthal mode number 
Pursuing further the comparison with the experimental results, we examine whet her the 



variou s spectral bumps correspond to specific azimuthal mode numbers. As in ^Schmitt et al. 



(|2008[ ), we correlate the signals computed at the same latitude (41°) but 72° apart in 
longitude. The signals are first narrow-band filtered, and we plot in figure [SJd the ampli- 
tude of the covariance (colour scale) as a function of the peak frequency of the filter, for 
time-delays between the two, converted into azimuthal mode number m (y-axis). As in 
the experiments (see figure [TJd), we find that each spectral bump corresponds to a sin- 
gle dominant (here negative by convention) azimuthal mode number m. The successive 
bumps have increasing azimuthal m (1, 2 and 3). 



3.3. Full Fourier transform 

In the numerical simulations, we can construct frequency spectra for each m. When the 
stationary regime is reached, we record 900 snapshots of the full fields, regularly spaced 
in time during 100 rotation periods: F(r, 6>, (/:?, t), where F can be either u or B. A two 
dimensional Fourier transform in the azimuthal and temporal directions ip and t gives us 
a collection of complex vectors representing the field for azimuthal number m and 
discrete frequency tx7, such that 

F(r,0,^,t) = ^^F-(r,e)e^(™^— (3.1) 
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Figure 7. Partial energy frequency spectra 8m (-uj) of the magnetic field at the outer surface 
(1 < r < 1.024, 0.35 < sin^ < 0.6) for different m. a) Equatorially symmetric part, b) Equato- 
rially- antisymmetric part. Frequencies are normalized by /, the rotation frequency of the inner 
sphere, and are shifted according to equation 13.31 Pm — 10~^, Re — 2 611 and A = 3.4 x 10~^. 



Note that the sign of the frequency has thus a precise meaning: positive (negative) 
frequencies correspond to prograde (retrograde) waves or modes. 
This allows us to compute partial energy spectra 

Sm{w) = r r ||F-(r,^)||' r smOdOdr. (3.2) 

Magnetic partial energy spectra fm(^) for the inner region (0.35 < r < 0.55, < 6 < 
tt) are shown in figure [6^ for m = to 4. They are dominated by a single peak, which 
moves towards positive (prograde) frequencies as m increases. This can be interpreted as 
the advection of stationary or low-frequency structures by a prograde fluid velocity. We 
therefore shift the frequency of the spectra in figure [6)3, according to: 



^shifted = — mf fluid (3.3) 

A good alignment of the spectral peak is obtained for f fluid — 0.5, which is compatible 
with the bulk fluid velocity in the outer region beneath the boundary layer (see figure 
[3^). This shift explains the linear evolution of the frequencies of the spectral bumps with 
m observed both in the DTS experiment (figure [TJd) and in the simulation (figure [SJd). 
It means that the peaks are caused by the advection of periodic structures by the mean 
flow, or by a non-dispersive wave. 

We now turn to the partial energy spectra of the magnetic field at high latitude (0.35 < 
sin^ < 0.6), at the surface of the sphere (1 < r < 1.024), displayed in figure [7^ andlTt 
for the equatorially symmetric and anti-symmetric parts, respectively. The frequencies 
are again shifted according to equation 13.31 This time, three peaks dominate the m = 
spectra. The spectrum is symmetric with respect to tu = since there cannot be prograde 
or retrograde propagation for m = 0: only latitudinal propagation or time-oscillations are 
permitted. The lateral peaks yield a frequency 1.8. As m increases, the lateral peak 

becomes dominant in the prograde direction, while it vanishes in the retrograde direction. 
We note that the peaks are well aligned in these shifted representation, meaning that 
these secondary fluctuations are also advected at roughly the same angular velocity as 
the central peak. But both stationary and propagating waves are required to explain 





-0.48-0.32-0.16 0.00 0.16 0.32 0.48 -0.036-0.024-0.012 0.000 0.012 0.024 0.036 



Figure 8. m = 2 mode structure given by the full Fourier transform method, for frequency 
uj'^ = 0.9, which corresponds to the m — 2 peak in figure [6^. Pm 10"^ Re 2 611 and 
A = 3.4 X 10~^. Energy densities are normalized by E%;. a) and b) isovalues in the cp — 0° 
meridional plane of the radial velocity and radial magnetic field, respectively, c) and d) mapview 
of the isovalues at r = 0.95 of the radial velocity and radial magnetic field, respectively. 



that this peak is not at zero frequency, and that it has both a prograde and a retrograde 
signature, and that the former dominates for m 7^ 0. 

Note that these secondary peaks do not show up in the regular frequency spectra or 
m-plots of point-measurements (figure [5]). This illustrates the interest and potential of 
the full Fourier transform method that we developed. 

3.4. Mode structure 

Picking the frequency that yields the maximum spectral energy density for a given m, 
we derive the structure of the corresponding mode. One example for m = 2 is shown in 
figure [HI where we plot, for both Ur and B^^ the structure of the mode in a meridional 
plane and in map view at r = 0.95. We selected a mode for which Br is symmetric with 
respect to the equator (and thus Ur is anti-symmetric). 

The meridional map for Ur reveals structures in the outer region, while Br shows 
similar patterns that extend deeper down to the inner sphere. While the map view of 
Ur near the outer boundary displays a rather chaotic structure, we find it remarkable 
that the structure of Br is very smooth and very similar to those retrieved in the DTS 



experiment, and well-accounted for in the linear modal approach of lSchmitt et al\ (|2Q12l ). 



4. Fluctuations and instabilities 

Having shown that our numerical simulations recover the essential features of the 
modes and spectra of the DTS experiment, even though their Reynolds number is much 
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Figure 9. Time- and- azimuth- averaged meridional map of the energy density of the fluctuations, 
a) Kinetic energy, b) Magnetic energy. Pm = 10""^, Re = 2 611 and A = 3.4 x 10~^. The energy 
densities are normalized by = |p (27r/)^ rf . 




Figure 10. Time- and-colatitude- averaged energy density of the fluctuations as a function of 
radius r. a) Kinetic energy, b) Magnetic energy. The energy densities are normalized by E^. 
The plots compare the fluctuations obtained for three simulations with the same Hartmann 
number (Ha — 16) and increasing Reynolds numbers {Re — 2 611, 5 070 and 10 100), and three 
simulations with the same Reynolds number {Re — 2 611) and increasing Hartmann numbers 
{Ha — 16, 19 and 31). Other dimensionless numbers as in Table [U Note that the magnetic 
energy is smaller than the kinetic energy by three orders of magnitude, but increases strongly 
toward the inner sphere for the simulation with the highest Reynolds number. 

smaller, we now examine where and how the modes are excited. The first guide we use 
is the location of the largest fluctuations. 

4.1. Energy fluctuations 
We compute the kinetic energy density as 5Ek = \p{{u — {u))'^) and the magnetic energy 
density as 5Em = 2^((^ ~ i^))'^)^ where () denotes time-averaging. We normalize both 
by a reference kinetic energy density E^^ = ^pl^^rf, and we integrate over azimuth. 

Figure [9] displays the resulting kinetic and magnetic energy densities of the fluctuations 
in a meridional plane for our reference simulation. We observe that the kinetic energy is 
maximum in the outer boundary layer, while the (much weaker) magnetic energy extends 
all the way to the inner sphere. 

Figure [To] shows the radial profiles obtained after integration over the colatitude 6. It 
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illustrates the effect of varying the Reynolds and the Hartmann numbers of the simula- 
tions. The fluctuations remain strongest in the outer boundary layer, but extend deeper 
inside the fluid with increasing Reynolds number. For the highest Reynolds number 
{Re = 10100), this causes the magnetic energy to strongly increase with depth, as the 
fluctuations interact with the larger imposed magnetic field near the inner sphere. 

4.2. Origin of the fluctuations 

It is beyond the scope of this paper to characterize the complete scenario by which 
instabilities develop in our geometry. However, we find it important to identify where 
the instabilities originate in order to understand the excitation of the modes we observe 
and extrapolate to other situations. The energy maps (figure [9j) as well as movies of the 
simulations (see online supplementary material) strongly suggest that fluctuations are 
initiated in the outer boundary layer. There is a large azimuthal velocity drop across the 
outer boundary layer, from the vigorously entrained fluid in the core flow to the outer 
container at rest. 

A detailed inspection of the numerical simulations reveals two types of instabilities, 
which do not occur in the same region but appear to be coupled through the meridional 
circulation: i) instability of a Bodewadt type layer at high latitude; ii) secondary insta- 
bility of a centripetal jet at the equator. We use the meridional snap-shots of figure [11] 
to illustrate these two mechanisms. 

4.2.1. Bodewadt layer instability 

The flow that appears when a fluid rotates at constant angular velocity above a flat 
disk at rest has been studied bv lBQdewadtl(|l940[ ) who has found the analytic expression 



of the boundary layer that develops at the surface of the disk. It is characterized by a 
large overshoot in t he azimuthal velo city profile caused by the centripetal radial circu- 
lation. As shown bv iLingwoodI ( 19971 ). this boundary layer is par ticularly un s table, and 



several teams ha ve analyzed the instabil ities that take pla ce fe.pi:.. ISayaslll987l : iLingwood 



19971 : iLopez et a l 2009; Ga uthier et a/.l ll999: Schouveiler et a/.ll2001f ). Two types of in- 



stabilities have been reported: axisymmetric rolls that propagate inwards (following the 
centripetal circulation of the boundary layer), and spiral rolls. A Bodewadt-type situa- 
tion is encountered in our geometry at high latitude. Figure [3^ shows a clear overshoot 
of the angular velocity at latitudes above about 40°, linked to a polewards meridional 
circulation. In this region, we observe polewards propagating axisymmetric rolls in our 
simulations when the inner sphere is spun from rest. This is best seen in the movies pro- 
vided as supplementary material online, but the signature of the rolls is clearly visible 
at high latitude in the three snap-shots of figure [Til 

It also shows up in the spatio-temporal representations of the instabilities at r = 0.95, 
just beneath the outer boundary layer (figure [12]) . Figure [T2k gives the axisymmetric 
part of ue as a function of time (x-axis) and latitude (^-axis) from — 7r/2 to 7r/2. The 
high latitude rolls show up as successive inclined lines in this time-latitude plot. 

Note that similar high latitude rolls have been reported in axisymmetric sim ulations 
of magnetized spherical Couette flow ( Hollerbach et a/.ll2007HBrito et a/.ll201l[ ). 



4.2.2. equatorial centripetal jet instability 

A different kind of instability takes place at the equator. In figure [T2k . polewards 
migrating rings yield a butterfly pattern, which also reveals that the equatorial down- 
welling instabilit y creates a meridi onal circulation of opposite sign around the equator 
for t < 30 turns. Brito et a/1 ( 201ll ) report this equatorial counter-rotating cell for some 



parameters in their axisymmetric equatorially- symmetric simulations. As shown in figure 
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Figure 11. Meridional snap-shots of the angular velocity of the fluid at successive times t (given 
in number of turns of the inner sphere) in our reference simulation. Pm — 10~^, Re — 2 611 and 
A = 3.4 X 10"^ a) t = 17.5; b) t = 93.8; c) t = 675. 



[TTk . it is associated to a sheet that draws fluid -and reduced angular momentum- from 
the outer boundary inside the sphere, in the eq uatorial plac e. It can probably be de- 
scr ibed as a c e ntrifu gal Taylor- Gortler instability ( Saricl[l993 ). similar to those observed 
by lNoir et all (|2QQ9[ ) in libration-induced flows in a sphere. As usually happens for these 
vortices, the non-linear evolution of the instability leads to mushroom- type downwellings 
(figures [TTbx and online movies (supplementary material)). 

At time t o:^ 24 turns, both the equatorial symmetry and the axisymmetry are broken 
by an m = 3 undulation, which rapidly disrupts the pattern of the fluctuations. Note 
however that axisymmetric bursts persist throughout, with amplitudes comparable to the 
initial ones. They propagate mostly polewards, but some occasionally cross the equator. 

The m = 3 undulation is best observed in the spatio-temporal plot of figure [T2b . 
which displays the non- axisymmetric fluctuations of the azimuthal velocity at a latitude 
of 10°, as a function of time and longitude (?/-axis) from to 27r. Until t 24 turns, 
there is no non- axisymmetric fluctuation, but dit t 24 turns an m = 3 mode appears 
(there are three maxima on a vertical line for a given t). After a few turns, this initial 
m = 3 undulation is replaced by chaotic fluctuations with dominant m = 1 and m = 
2 contributions, which travel in the prograde direction with approximately the same 
velocity (given by the slope of the color streaks in this figure). 

The m = 3 secondary in s tability is similar to those observed in non-magnetic spher- 
ical Couette ( iDumasI Il99ll : iGuervillv fc CardinI [ioTol ) or with an axial magnetic field 
(|Hollerbaclj |2QQ9[ ). There, it takes place on the centrifugal equatorial jet, which is a 
primary feature of these flows. 

In our case note that, while the equatorial counter-rotating cell is essential for the 
centripetal jet to form, the time-averaged meridional circulation (shown in figure [Sb) 
does not show this feature, as if the interplay of the developed instabilities had erased it. 

4.2.3. threshold of instability 

Although we do not intend to decipher the complete scenario of instability, we have 
determined the threshold of instability, which is found at Rcc = 1 860, with a critical 
azimuthal mode number of 2. Interestingly, it looks that the two (coupled) instabilities 
described above are present from this threshold. We can relate this threshold to the 
critical Reynolds number of the boundary layer. Following iLingwoo 3 (|l997h . we define 
the local Reynolds number re = uo'^s^l/v^ where / is the thickness of the laminar boundary 
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Figure 12. Spatio-temporal representation of the velocity fluctuations in the simulation at 
radius r = 0.95. Pm = 10~^, Re — 2 611 and A = 3.4 x 10~^. a) ug isovalues of the axisymmetric 
flow (m = 0) in a latitude-time plot, b) isovalues at 10° latitude in a longitude-time plot 
(m = term removed) . Velocities are normalized by the reference velocity Qri . Time is given in 
rotation periods. 



layer: / = ^Jv|bf^ ^ with b:^ the dimensional angular velocity of the fluid with respect to 
the wah at a position specified by its dimensional cylindrical radius . We can relate it 
to our global Reynolds number Re by: 



re = (jjs 




(4.1) 



where the angular velocity cj just outside the boundary layer is adimensioned by 2i\j and 
the cylindrical radius s by ro, as before. Picking u 0.7 at a latitude of 45° (s 0.7) 
from fig ure [H we get re^ 36. This is somewhat larger than the critical value of 21.6 
found bv lLingwoodI ( 19971 ) for the absolute instability of a pure Bodewadt layer, for which 
she predicts a critical mode number rric = ^c^^c = —0.1174 x 21.6 —2.5, which is not 
incompatible with our observation of an initial m = 2 or m = 3 pattern. It is difficult to 
assess wh ether the higher thresho ld we get is due to a stabilizing effect of the magnetic 
field as in oresco fc Alboussierd (!M04), or to the spherical geometry. 



In any case, all the experiments analyzed bv lSchmitt et al\ (|2008L l2012f ) are far above 
this threshold. Our main conclusion at this stage is that the fluctuations we observe 
initiate in the outer boundary layer, where the influence of the magnetic field is probably 
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Figure 13. Energy of the fluctuations in the DTS experiment (all energy densities are normal- 
ized by E'^ = ^pQ?r1). a) Selected radial profiles of the kinetic energy of the fluctuations 5Ek 
deduced from ultrasonic Doppler radial velocity measurements at latitude —20°. b) Selected 
radial profiles of the magnetic energy of the fluctuations 5Em deduced from measurements in 
a sleeve at two different latitudes (10° -solid green- and 40° -dashed blue). The symbols corre- 
spond to different rotation rates / of the inner sphere, given in Hz in the legends, c) Evolution of 
the kinetic ('x') and magnetic ('+') energies as a function of the inverse of the Elsasser number. 
The black symbols are from the numerical simulations. The colored symbols are deduced from 
the maximum value of the experimental radial profiles of 5Ek and 5Em at various latitudes and 
forcings /, using the same markers. Note that magnetic energies are multiplied by 100. 



negligible. Because the fluid is in rapid rotation beneath the outer sphere at rest, the outer 
boundary is very unstable, and subject to non-geostrophic instabilities. We therefore 
expect a radically different behaviour when the outer sphere spins and the boundary 
layer is of Ekman type. 



4.3. Comparison with experimental results 

We cannot measure the total kinetic and magnetic energies of the fluctuations in the 
DTS experiment. However, we can get a quantitative assessment of the energy of the 
fluctuations as a function of radius, at given latitudes. The kinetic energy is obtained 
from the fluctuations of the radial velocity measured by ultrasound Doppler velocimetry 
along a radial shot at a latitude of —20°. The magnetic energy is derived from the 
fluctuations of the azimuthal component of the magnetic field measured at three latitudes 
10°, 20° and 40° and at 6 different radii, using Hall probes inserted in a sleeve, after 
removing a contribution at the rotation frequency / and harmonics, which is due to small 
heterogeneities of the imposed magnetic field. All energy densities are scaled with = 
^pVt^rl. As for the simulations, we integrate over azimuth by multiplying the measured 
rms by 27rrsin6> and convert to energy. In order to relate to the numerical results, we 
assume that fluctuations are isotropic. Additional measurements of the magnetic energy 
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Figure 14. Same as figure[5]but for a simulation where the Lorentz force is nulled out for m / 0. 
Sharp spectral peaks are observed on all three components of the magnetic field. Neighbouring 
frequencies can have different azimuthal modenumber m. 



from radial and orthoradial probes partly support this hypothesis. Nevertheless, the 
comparison remains approximative. 

The kinetic energy profiles (figure [T3k ) confirm that fluctuations are strongest near 
the outer surface. The maximum is deeper than in the numerical simulations (compare 
with figure [TOk. a consequence of the much higher Reynolds number. Note however that 
the thin viscous boundary layer cannot be resolved from the Doppler velocity profiles. 

The magnetic energy profiles (figure [T3b) clearly show that fluctuations are strongest 
near the inner sphere. Figure [TOb shows that only the simulation with the highest 
Reynolds number displays this behaviour. 

Figure [T3h compares the kinetic and magnetic energies of the fluctuations obtained 
from both the simulations (black crosses) and the experiments (colored symbols). Since 
we don't have the full latitudinal dependence in the experiments, we simply take the 
maximum of each profile as an estimate of the overall energy. The horizontal axis is 
A~^, the inverse of the Elsasser number. It measures the ratio of the inertial force to 
the Lorentz force. The magnetic energy remains much smaller than the kinetic energy. 
It clearly increases with in the simulations: velocity fluctuations penetrate deeper 
into the fluid and induce larger magnetic fluctuations because the imposed magnetic field 
is stronger there. The experimental data follow the same trend for small forcing /, but 
there seems to be a strong drop near = 100, before it increases again. We have no 
explanation for this behaviour. 

4.4. The role of the Lorentz force 

Although magnetic energies are much smaller than kinetic energies, the Lorentz force 
plays a major role. The strong imposed dipolar magnetic field governs the dynamics of 
the mean flow in the DTS experiment. In particular, the very efficient entrainment of 
the fluid by the conductive inner sphere, and the zone of super-rotation next to it, are 
entirely due to the presence of the magnetic field. 

In order to see the effect of the Lorentz force on the fluctuations, we have run a 
simulation in which the Lorentz force is nulled out except for m = 0. We find that the 
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Figure 15. Time- and-colatitude- averaged energy density of the fluctuations as a function of 
radius r. This plot compares the kinetic (solid lines) and magnetic (dashed lines) energies of the 
reference simulation (in red) with a simulation (in blue) with the same parameters {Pm — 10~^, 
Re — 2 611 and A = 3.4 x 10~^), but in which the Lorentz force has been nulled out for m 7^ 0. 
Note that the magnetic energy is two orders of magnitude larger in the latter case. 



radial profile of angular velocity at the equator remains essentially the same, illustrating 
that non-linear interactions of fluctuations with m 7^ barely contribute to the mean 
flow. Frequency spectra of the surface magnetic field (figure [T4k ) still display spectral 
bumps, but they are narrower and much more intense. Furthermore, the azimuthal mode 
number analysis (figure [Mb) reveals that modes of a given m show up at several distinct 
frequencies. 

The damping effect of the Lorentz force is best illustrated by plots of the radial profile of 
the energy of the fluctuations. Figure [15] displays the resulting radial profiles of the kinetic 
and magnetic energies, integrated over time and colatitude. While the kinetic energy of 
the fluctuations is negligible for r < 0.6 when the Lorentz force is present, fluctuations 
invade the complete fluid shell when it is suppressed. Fluctuations are largest beneath 
the outer shell where the magnetic field is weakest, but even there fluctuations are much 
weaker when the Lorentz force is active. 

The magnetic energy of the fluctuations reaches only a thousandth of the kinetic 
energy in the DTS configuration. Not surprisingly, when the Lorentz force is suppressed, 
it jumps by a factor of 100 near the inner sphere, where the imposed magnetic field 
is strongest. This demonstrates how dangerous it can be to infer magnetic energy and 
dissipation from flo w solutions compu ted without th e feedb ack of the Lorentz force (see 
discussion between Glatzmaier (2008) and Liu et al. ( 20081 )). 



5. Discussion 

We have obtained bumpy frequency spectra in numerical simulations of the magnetized 
spherical Couette flow (figure [5^) . They compare very well with spectra obtaine d in the 
DTS experiment from magnetic and electric time-records (^Schmitt et "^ l2008l) . As in 
the experiment, the dominant azimuthal mode number is incremented by 1 from one 
frequency bump to the next (figure [Sb). Very large Reynolds numbers are not required 
for getting this behaviour, but one needs to accumulate long time-series (typically 300 
turns of the inner sphere) for the peaks to show up clearly in the spectra. However, we 
note that bumps in the numerical spectra are not as pronounced as in the experiments. 

We hav e developed a new method, which bridges the gap between the line ar modal 
approach ( Rieutord fc Valdettarol I1997I : iKellev fi^ a/.ll2007l : ISchmitt et fl/.ll2008l ) and full 
non-linear simulations and experiments. By performing a time-domain Fourier trans- 
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form of the full fields for each azimuthal number m, we recover the dominant frequencies 
(figures [6] and [7j) and obtain the structure of the modes (figure [8]), which can then be com- 
pared to Hnear solutions and to experimental observations. We think that t his approach 



will h elp identifying the mode selection mechanism in other experiments (|Kellev et al 

hold ). 

Snap-shots (figure [TT]) and spatio-temporal plots (figure [T2]) reveal a rather different 
story, in which chaotic instabilities are swept by the fiow. We could show that these two 
views are dual: since the instabilities circle around the sphere, the parts that are in phase 
between two successive passages are statistically enhanced. Since the spectral bumps are 
more pronounced in the experiments, this effect appears to be more efficient at large 
Reynolds number. 

The maps of the kinetic energy of the ffuctuations (figures [9^1, flOk and [T3k ) indicate 
that they initiate in the outer boundary layer, with only minor infiuence of the magnetic 
field. Instabilities appear above a critical Reynolds number Rec = 1 860. We identify two 
types of instabilities: i) axisymmetric polewards migrating rolls at high latitude; ii) non- 
axisymmetric (m = 2 at the threshold) secondary instabilities of an equatorial centripetal 
jet. The first type is similar to the instabilities of a Bodewadt layer. The second type 
resembles the jet-instability of centrifugal equatorial sheet in non-magnetized spherical 
Couette fiow. The two instabilities are coupled by the meridional circulation (they trigger 
one another), and the system quickly evolves towards a chaotic state in which outer 
boundary layer instabilities are swept around by the azimuthal and meridional large 
sc ale fiows. 

ISchmitt et al\ ( 2008 ) observe that fiuctuations in DTS are delayed and reduced when 
the outer sphere is spinning. We think that this is bec ause the bounda ry layer is then 
closer to an Ekman-type, which is much more stable (|Lingwoodl Il997l ) . and that non- 
geostrophic instabilities are hampered. 

The fiuctuations of kinetic energy are much larger than that of magnetic energy (figure 
[T3h). which are mostly slave of the former. If we assume that the magnetic fiuctuations 
scale 8iS b aRmB, we obtain that the ratio of the magnetic to kinetic energy behaves 
as SEm/^Ek — o?Lu^^ where Lu is the Lundquist number defined in section [2^ Since 
the Lundquist number is of order 1 in both the experiments and the simulations the a 
pre-factor must be rather small to explain the observed energy contrast. As a matter of 
fact, direct measurements of the mean induced azimuthal magnetic field yield a o:^ 0.1. 
Both the DTS measurements (figure [T3b) and our largest Reynolds number numerical 
simulation (figure [TOb) display a strong increase of the magnetic energy fiuctuations when 
getting closer to the inner sphere. This appears to be essentially the consequence of the 
strong increase of the imposed dipolar field there. Indeed, the local Lundquist number 
increases from Lii = 0.5 at the equator of the outer sphere to Liz = 12 at the equator of 
the inner sphere. 

At first order, we expect both energies to be proportional to the square of the im- 
posed inner sphere velocity. However, we note that when scaled accordingly, the kinetic 
energy tends to decrease when the forcings; is increased, while the scaled magnetic en- 
ergy increases (figure [T3b). iBrito et o/l ( 2011 ) showed that the energy of the mean fiow 
behaves similarly, and proposed that this is a consequence of the increasing turbulent 
friction at the outer surface: as the friction increases, the core-fiow is slowed down, while 
the shear between the spinning inner sphere and the fiuid increases, inducing a stronger 
magnetic field. We think that another effect explains the trend observed for the energy 
fluctuations: as the forcing increases, the damping effect of the magnetic fleld decreases. 
Instabilities penetrate deeper into the fluid and produce larger magnetic fluctuations, 
even though their scaled kinetic energy is reduced because of the decreased velocity drop 
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across the outer boundary layer. We note that for large forcing / the magnetic energy 
is smaller in the experiments than in the simulations, suggesting again the role of the 
strong turbulence. 

Although the magnetic energy is very small, the Lorentz force plays a major role: 
it determines the very efficient entrainment of the fluid by the spinning inner sphere, 
but it also heavily damps the fluctuations in most of the fluid. When we remove the 
Lorentz force for m 7^ 0, fluctuations invade the fluid (figure [15]) , and sharper and more 
numerous frequenc y peaks are o bserv ed in the spectra (figure [T4|) . This gets closer to 
the observations of lKellev et al\ ((2007D, where the imposed magnetic field was weak and 



only served as a n iarker of the flow. It could also explain the sharp W modes reported 
bv lSchmitt et all ( 2008 ) for the largest rotation rates in the DTS experiment. 



One could have expected instabilities to develop in the inner region near the equator, 
where the flow obeys Ferraro law, and where a small velocity perturbation produces a 
large Lorentz force. This does not appear to be the case. In the present study, we have 
kept the Lundquist number small, as in the DTS experiment. Alfven waves are therefore 
damped out rapidly. They might still contribute to shaping the modes near the inner 
sphere. It would be interesting to investigate the turbulent regime in the DTS geometry 
at larger Lundquist number. 

We thank D. Jault for stimulating discussions. R. Hollerbach and three other refer- 
ees helped us improve our manuscript. We gratefully acknowledge the support of CNRS 
and Universite de Grenoble through the collaborative program "Turbulence, Magneto- 
hydrodynamics and Dynamo" . Part of the numerical simulations were run at the Service 
Commun de Calcul Intensif de lObservatoire de Grenoble (SCCI). 

Supplementary movies are available at journals.cambridge.org/flm. 

Movie 1. Map view of the time evolution of the azimuthal velocity beneath the surface of 
the outer sphere (r = 0.95) in the reference numerical simulation (Pm = 10~^, Re = 2 611 
and A = 3.4 x 10~^). At the origin time, the fluid is at rest and the rotation rate of the 
inner sphere is set to /. The frame index is displayed. There are 6 frames per turn. The 
movie lasts 100 turns. Note that the first instabilities appear at the equator and are 
axisymmetric. After about 24 turns, non-axisymmetric instabilities show up. 

Movie 2. Time evolution of the angular velocity in a meridional plane {(f = 0) for the 
same simulation as in movie 1. 
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